On the use of Kolmogorov-Landau approach in deriving various correlation functions 

in 2-D incompressible turbulence 
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We look at various correlation functions, which include those that involve both the velocity and 
the vorticity fields, in 2-D isotropic homogeneous decaying turbulence. We adopt the more intuitive 
approach due to Kolmogorov (and subsequently. Landau in his text on fluid dynamics) and show 
that how the 2-D turbulence results, obtainable using other methods, may be established in a simpler 
way. Also, some experimentally verifiable correlation functions in the dissipation range have been 
derived for the same system. The paper also showcases the inability of the Kolmogorov-Landau 
approach to get the "one-eighth law" in the enstrophy cascade region. As discussed in the paper, 
this may raise the spectre of logarithmic corrections once again in 2-D turbulence. 

PACS numbers: 47.27.i,47.10.ad 



I. INTRODUCTION 

Kolmogorov's four-fifths law[l] is a landmark in the 
theory of turbulence because it is an exact result. In 
three spatial dimensions, this law says that the third or- 
der velocity correlation function behaves as: 



n 3\ 



{v{x + r) — v{r)} . 



where e is the rate per unit mass at which energy is 
being transferred through the inertial range. The iner- 
tial range is the intermediate spatial region postulated 
by Kolmogorov where the large scale disturbances (flow 
maintaining mechanisms) and the molecular scale viscous 
dissipation play no part. This result is of such central sig- 
nificance that attempts are regularly made to understand 
it afresh and to extend it in other situations involving 
turbulence. There appears to be following two primary 
methods of obtaining this result: 

1. The original Kolmogorov method put forward in 
details in the fluid dynamics text due to Landau 
and Lifshitz[2]. There is no external forcing in this 
approach and the equality of dissipation rate and 
forcing rate for the energy is never enforced. 

2. A "field-theoretic technique" because of its analogy 
to anomalies in the calculations of certain field- 
theoretic correlation functions. This is the path 
following in the more recent turbulence text by 
Frisch[3]. In this approach, there is an external 
forcing that maintains a steady state turbulence. 

The two approaches yield the same result as they should. 
If we consider the two-dimensional turbulence, then in 



the inviscid limit, we have two conserved quantities - en- 
ergy and enstrophy. This gives rise to two fiuxes with the 
enstrophy flux occurring from the larger to the smaller 
spatial scales. The energy flux goes in the reverse direc- 
tion. Recently, Bernard [4] has used the second of the 
above mentioned techniques to obtain the third order 
structure function for both the energy and the enstrophy 
cascade regions. Interestingly enough, he finds an exact 
answer for the two-point third order correlation function 
in the enstrophy cascade. Now, the correlation functions 
in the enstrophy cascade are often plagued by logarith- 
mic corrections [5] and in view of this Bernard's result 
is particularly interesting. We believe that the issue is 
important enough that a derivation of Bernard's result 
using the Kolmogorov-Landau approach should be use- 
ful. To the best of our knowledge and literature survey, 
such an approach remains unreported. So this is what we 
have attempted here and while the energy cascade result 
does come out in agreement with Bernard's, it appears 
that the Landau treatment cannot yield the prefactor of 
the so-called "l/8th-law" for the two-point third order 
correlation function in the enstrophy cascade dominated 
regime. This would be the case if logarithmic corrections 
exist. We believe that this raises the spectre of logarith- 
mic corrections once again. Besides, we also have derived 
some other correlation functions which deal with vorticity 
fields in the inertial region and also some two-point sec- 
ond order correlation functions in the dissipative region 
following the arguments of Landau, thereby consolidating 
the equivalence between the two approaches mentioned 
in the beginning. 



II. SECOND ORDER CORRELATION 
FUNCTION FOR ENSTROPHY CASCADE IN 
DISSIPATION RANGE 
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It is a well-established fact that there exists a direct- 
cascade of enstrophy in two-dimensional (2-D) turbu- 



lence. One defines enstrophy as F = i J, 



,2^2; 



2 J all space 
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where uj = dj.Vy — dyV^ is the vorticity in the Carte- Again, having assumed incompressibility, one may write: 



sian coordinates; v being the velocity field. As we shall 
consider incompressible fluids only (V.tT = 0), we shall 
take density to be unity and let p take over the task of 
representing position vector in 2-D plane. The enstro- 
phy flows through the inertial range and gets dissipated 
near dissipation scale. Using the antisymmetric symbol 
EajS that has four components, viz. = £22 = and 
£12 = —£21 = 1, one may define the mean rate of dissi- 
pation of enstrophy per unit mass as: 



T] = i/(Va;.Vw) 

ry = veTa£e0{{drd^Va){d0d^V(})) 



(1) 



Here, angular brackets denote an averaging procedure 
which averages over all possible positions of points 1 and 
2 at a given instant of time and a given separation. Now, 
if Vi and V2 represent the fluid velocities at the two neigh- 
bouring points at pi and p2 respectively, one may define 
rank two correlation tensor: 



BaP = {{via - Vla){V2l3 - Vif))) 



(2) 



For simplicity, we shall take a rather idealised situation 
of turbulence fiow which is homogeneous and isotropic on 
every scale, a case achievable in practice in a vigorously- 
shaken-fluid left to itself. The component of the corre- 
lation tensor will obviously, then, be dependent on time, 
a fact which won't be shown explicitly in what follows. 
As the features of local turbulence is independent of av- 
eraged flow, the result derived below is applicable also to 
the local turbulence at distance p much smaller than the 
fundamental scale. Isotropy and homogeneity suggests 
following general form for i3a/3[6] 



B^p = Ai{p)Sc,0 + A2{p)p°„t 



(3) 



where Ai and A2 are functions of time and p. The Greek 
subscripts can take two values p and _L which respectively 
mean the component along the radial vector p and the 
component in the transverse direction. Einstein's sum- 
mation convention will be used extensively. Also, 

P = P2-P1, p°^ = Pa/\p\, P°p = i, p1=0 

using which in the relation (3), one gets: 

= - P>^) + S^,P>^ (4) 

One may break the relation (2) as 

Ba0 = {viaVl/}) + {V2aV20) - {viaV2p) - {v2aVl0) (5) 

and defining 

baf3 = {viaV2f3) (6) 

one may proceed, keeping in mind the isotropy and the 
homogeneity, to write 



Bal3 = {v'^)Saf} - '^bal3 



(7) 



dpBa/s = 



B',, + -{B,,-B^^) = 



B 



pB'pp + B, 



pp 



(8) 



where the equation (4) has been used and prime (') de- 
notes derivative w.r.t. p. Near the dissipation region the 

flow is regular and its velocity varies smoothly which al- 
lows to expand w in a series of power of p. One must take 
V ^ p^ neglecting the higher powers ~ p is not taken 
because it leads to the contradictory result that = 
as can be seen from the relation (1)). So, treating o as 
a proportionality constant, let Bpi, = ap'^, which means 
B±± = 5ap'^ (using equation (8)) and hence, 

1 5 

{vicxV20) = -^{v'^)Saf3 - -^ap^Safi + 2ap'^ Pap0 

^ {{diTdi^via){d2ed2yV20)) = -72a50r5a0 + 24:a500dar 

+ 24:aSa05/}r 

^ £TaSei3{{dTd^Va){d0d^V0)) = -192a (9) 



Bnn = — 



VP 

192z/ 



(10) 



In the equation (10), we have put pi « p2, for these 
relations are assumed to be valid for arbitrarily small p. 
While writing the relation (10), equation (1) has been re- 
called. This Bpp is the second order correlation function 
for enstrophy cascade in dissipation range. 



III. THIRD ORDER CORRELATION 
FUNCTION FOR ENSTROPHY CASCADE IN 
INERTIAL RANGE 

Let's again define: 

&a/3,7 = {viaVl/3V2'y) 

Invoking homogeneity and isotropy once again along with 
the symmetry in the first pair of indices, one may write 
the most general form of the third rank Cartesian tensor 
for this case as 

&a/3,7 = C{p)SapP°+D{p){S^i3p°^ + SajP°^) 

+F{p)plp'>pp° (11) 

where, C, D and F are functions of p. Yet again, incom- 
pressibility dictates: 

— Oa/3,7 = -Q-Oapa = 



d- 



27 



C, 2D ^ 2D' 2D 

Oal3 H Oa/3 H Oa/3 H n-paPp o-PaPjS 

p p p^ P 

F' F 

+ ^PaPl3 + -^PaP0 = 0(12) 
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Putting a = (3 in equation (12) one gets: 

^ „ „ constant 
2C + 2D + F = = 



(13) 



where, it as been imposed that bap^j should remain finite 
for p = 0. Again, using equation (12), putting a ^ (i and 
manipulating a bit one gets: 



D = -\{pC' + C) 



(14) 



using which in relation (13), one arrives at the following 
expression for F: 



F = pC' -C 



(15) 



Defining 



= 2(6a/3,7 + ^7/3,a + ^a7,/3) (16) 

and putting relations (14) and (15) in the equation (16) 
and using relation (11), one gets: 

-Ba/37 = -^pO'{Sa0p" + S^0p^ + Sa-yP^p) 

+ &{pC' -C)pIp''0P° (17) 
-6C (18) 



^ 5, = B, 



ppp 



which along with relations (14), (15) and (11) yields the 
following expression: 



S3 1 

&c</3,7 = ^Sapp" + Y^ipS's + S3){S^0Pa + ^aj, 

-J(P^3-'53)P>K 

Navier-Stokes equation suggests: 

dt 
d 



(19) 



-Via = -Vl^di^Via - diaPi + vdi^.di^.Vi^ (20) 



dt 



V2f3 = -V2id2^V2j3 - ^2/3^2 + vd2^d2^V2(3 (21) 



multiplying equations (20) and (21) with V2(} and Via 
respectively and adding subsequently, one gets the fol 
lowing: 

d 



Using equations (7) and (19), one can rewrite equation 
(23) as: 



d ( as. 



p dp V dp 



(24) 



As wc arc interested in the cnstrophy cascade, the first 
term in the R.H.S. is zero due to homogeneity and the 
first term in the L.H.S. is zero because of energy remains 
conserved in 2-D turbulence in the inviscid limit (and it is 
the high Reynolds number regime that we are interested 
in); it cannot be dissipated at smaller scales. Also, as we 
are interested in the forward cascade which is dominated 
by enstrophy cascade, on the dimensional grounds in the 
inertial region Bpp (as it may depend only on j] and p) 
may be written as: 



(25) 



where ^ is a numerical proportionality constant. Hence, 
using the relation (25), the equation (24) reduces to the 
following differential equation: 



Qp^dp^'' ""''> ~ pdpY dp 



A 



(26) 



which when solved using relation(18) in the limit of infi- 
nite Reynolds number {v ^ 0), one gets 



Bppp = --^p 



(27) 



-dia{piV20) - d2(3{p2Vla) + vdi^di^{via 

+ vd2^d2^{viaV2p) 

Due to isotropy, the correlation function for the pres- 
sure and velocity {{P1V2)) should have the form f{p)p/\p\. 
But since, da{piV2a) = due to solcnoidal velocity field 
f{p)p/\p\ must have the form constant x (p/|p|^) that in 
turn must vanish to keep correlation functions finite even 
at p = 0. Thus, equation (22) can be written as: 



The relation (27) is the expression for the third order 
correlation function of the isotropic and homogeneous 2- 
D decaying turbulence in the inertial range of the for- 
ward cascade. The fact that wc could not determine 
an exact value for A is what is disturbing. Bernard's 
procedure leads to'A'=— 1/4. One knows that in two- 
dimensional turbulence the energy spectrum in enstrophy 
cascade regime is of the form Kr]'^^^k~'^{ln{k/ki)}^^/'^ 
(where K isa numerical constant and A; 3> fci) and hence, 
recalling that the integral of energy spectrum over k is 
directly related to Bap. in the relation (25), the dimen- 
sional form of Bpp should have a logarithmic correction. 
This would immediately lead to a logarithmic correction 
in Bppp, making the existence of 'A' as a constant mean- 
^^^^ingless. It is possible that the difference between the 
Bernard's and our approaches shows a basic difference 
^ 'between a forced turbulence and a decaying turbiilcnce 
which comes to the fore in the example being dealt with. 



d_ 
di 



bap = dj{baj^p + bp^^a) + 2vd^d^bap 



(23) 



IV. THIRD ORDER CORRELATION 
FUNCTION FOR ENERGY CASCADE IN 
INERTIAL RANGE 

Suppose in the homogeneous isotropic fully-developed 
turbulence in 2-D, energy is being supplied and the mean 
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rate of injection of energy per unit mass is denoted by e. 
Let us concentrate on the inverse energy cascade. Then 
technicahy we have to proceed as in the previous sec- 
tion and on doing so one would arrive at the differential 
equation (24); only that now the arguments would dif- 
fer. In the larger scales viscosity is not as significant and 
anyway we shall be interested in the infinite Reynolds 
number case which would mean that the last term in the 
R.H.S. of equation (24) would go to zero. One obviously 
would also set \^{v'^) = e and lets assume -^^Spp « in 
the inverse cascade regime, justification of which can be 
sought from the fact that the ultimate result that is ob- 
tained has been experimentally and numerically verified. 
So we are left with the following diff'erential equation: 



_i_d_ 
6p^ dp 



= £ 



Bppp = +—ep 



(28) 
(29) 



where in the last step the integration constant has been 
set to zero to prevent Bppp from blowing up at p = 0. 



V. SECOND ORDER VORTICITY 
CORRELATION FUNCTION IN DISSIPATION 
RANGE OF ENSTROPHY CASCADE 

Let us define: 

W = {1JJ1UJ2) (30) 
and, O = {{ijJ2 — iO\){ijJ2 — wi)) (31) 

Due to homogeneity, Q. may be expressed as: 

9. = 2{uj'^) - 2W (32) 

In the dissipation range: w ~ so, w ~ p and hence we 
may, choosing a proportionality constant b (say), assume: 



Using relations (30), (32) and (33), one gets: 



(33) 

(34) 
(35) 



=> {{di^Uj){d2c,ij)) = 2b 

But we know, 

ri = v{(dauj){dauj)) (36) 

So, relation (35) would yield: 

r? = 2ub (37) 

where, we have put pi « p2 in the relation (35), for, being 
in the dissipation range, these relations are assumed to 
be valid for arbitrarily small p. Using relations (33) and 
(37), one arrives at a experimentally verifiable result for 
two-point second order vorticity correlation function: 



(38) 



VI. THIRD ORDER MIXED CORRELATION 
FUNCTION IN INERTIAL RANGE OF 
ENSTROPHY CASCADE 

We start by defining a two-point third order mixed 
correlation tensor in inertial range: 

= {{V2l3 - Vli3){uJ2 - wi)(a;2 - Wl)) (39) 
^ Ct/3 = 2M/3 + AW/s (40) 



where. 



and, 



Wff = {V113UJ1LJ2) 

= {uJlUJlV2p) 



(41) 
(42) 



Let due to isotropy and homogeneity, we can write fol- 
lowing form for Mg: 



M0 

_d_ 

P2/3 

d_ 
dp 

M{p) 



= M{p)pl (43) 
Mfj = {u)iu)id2i3V20) = (44) 
M(p) 



M{p) 



= 



P 

constant 







(45) 
(46) 



In the relation (44), we are assuming incompressibility 

and in writing the relation (46) we have taken into ac- 
count the fact that Mfj should remain finite when p = 0. 
Relations (43) and (46) imply that: 



=0 

using which in the relation (40), we get: 



(47) 



(48) 



Prom the equations (20) and (21), we may write respec- 
tively: 



d_ 

di 
d_ 

di 



UJl — — Wl-y(?i-yCJl + vdl-ydl-yUJl (49) 

W2 = -l'27t?27W2 + vd2'yd2'yUJ2 (50) 



Multiplying equations (49) and (50) by uj2 and loi re- 
spectively and adding subsequently, we get the following 
differential equation: 



d_ 

dt' 



W = 2d0Wf} + 2vdBd3W 



(51) 



where we have used the fact dfs = —dip = 92/3- Using re- 
lations (32) and (48) in the equation (51), one gets for the 
inertial range for the enstrophy cascade in homogeneous, 
isotropic and fully-developed freely decaying turbulence 
in 2-D in the infinite Reynolds number limit {i.e., v ^ Q) 
following differential equation: 



dV ' 2dt 2p5p^^ 



-2r)p 



(52) 
(53) 
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In getting relation (53) from the equation (52). we have 
used the facts: ^^(f^^) = — ?? and ^-g^^ « as it may 
be supposed that the value of O varies considerably with 
time only over an interval corresponding to the funda- 
mental scale of turbulence and in relation to local tur- 
bulence the unperturbed flow may be regarded as steady 
which mean that for local turbulence one can afford to 
neglect in comparison with the enstrophy dissipation 
rate rj. This result (relation (53)) has gained importance 
by serving as the starting point in deriving various rig- 
orous inequalities for short-distance scaling exponents in 
2-D incompressible turbulence [7]. 

VII. CONCLUSION 

The Kolmogorov-Landau approach has been invoked 
in 2-D homogeneous isotropic decaying fluid turbu- 
lence to arrive at the various correlation functions 
earlier obtained using different methods. Also, some 



experimentally verifiable correlation functions in the 
dissipation range have been derived, which to the best 
of our knowledge, have not been reported anywhere 
in the literature yet. The results derived here are 
'exact' (though not rigorous) something which is a far 
cry in the literature on turbulence. Besides, we have 
shown that the approach we have adopted here possibly 
highlights the fact that the logarithmic corrections in 
2-D turbulence may hinder the proof of "one-eighth law" 
in the enstrophy cascade region, and thereby showcases 
the handicap of the Kolmogorov-Landau approach in 
2-D fluid turbulence. 
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